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ABSTRACT 

This paper presents deep high quality photometry of globular cluster (GC) systems belong- 
ing to five early-type galaxies covering a range of mass and environment. Photometric data 
were obtained with the Gemini North and Gemini South telescopes in the filter passbands 
g', r', and The combination of these filters with good seeing conditions allows an excel- 
lent separation between GC candidates and unresolved field objects. In fact, our previously 
published spectroscopic data indicates a contamination level of only ~10 percent in our sam- 
ple of GC candidates. Bimodal GC colour distributions are found in all five galaxies. Most 
of the GC systems appear bimodal even in the (g'-r') vs (r'-i') plane. A population of re- 
solved/marginally resolved GC and Ultra Compact Dwarf candidates was found in all the 
galaxies. A search for the so-called "blue tilt" in the colour-magnitude diagrams reveals that 
NGC 4649 clearly shows that phenomenon although no conclusive evidence was found for 
the other galaxies in the sample. This "blue tilt" translates into a mass-metallicity relation 
given by Z oc 7\/0-28±o.03 -pj^jg dependence was found using a new empirical (g'-i') vs [Z/H] 
relation which relies on an homogeneous sample of GC colours and metallicities. This paper 
also explores the radial trends in both colour and surface density for the blue (metal-poor) and 
red (metal-rich) GC subpopulations. As usual, the red GCs show a steeper radial distribution 
than the blue ones. Evidence of galactocentric colour gradients is found in some of the GC 
systems, being more significant for the two SO galaxies in the sample. Red GC subpopulations 
show similar colours and gradients to the galaxy halo stars in their inner region. A GC mean 
colour-galaxy luminosity relation, consistent with [Z/H] oc L gO-26±o.08^ jg present for the red 
GCs. An estimate of the total GC populations and specific frequency Sat values is presented 
for NGC 3115, NGC 3923 and NGC 4649. 

Key words: galaxies: elliptical and lenticular, cD -galaxies: star clusters: general -globular 
clusters: general 



1 INTRODUCTION 

The idea that globular clusters (GCs) are good tracers of very early 
events in the lives of galaxies was put forward many years ago 
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Council (Canada), CONICYT (Chile), the Australian Research Council f E-mail:favio@fcaglp.unlp.edu.ar 
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(e.g. lEggen, Lvnden Bell & Sandagelll962l : ISearle & Zinnlll978h . 
A thorough description of all the im pori:ant aspects is given, for 
example, in ISrodie^ Stradeil ll2006l) . However, some key issues, 
such as the relation between GCs and galaxy fie ld stars, still re- 
main as open questions (see iKissler-Pati j|200^ . If such a rela- 
tion does exist, some galaxy properties could be described in terms 
of the characteristics of their globu lar cluster systems (GCS). A 
tentative approach in this direction Jporte. Faifer & Geisleill2007l : 
iForte. Vega & Faifeill2009l) assumes that GCs formation efficiency 
depends on chemical abundance, and leads to a quantitative con- 
nection between field stars and GCs. In turn, this type of con- 
nection allows further analysis of the diffe rent GCs-galaxy forma- 
tion scenarios found in the literature ( e.g. Ashman & Zepf '1992'; 
iForbes et al.lll997l : iBeaslev et all2002l : iMuratov & Gnedin,201Q) . 

In this work we present new photometry for the GCS of five 
early-type galaxies, using the imaging capabilities of the GMOS 
instruments on the Gemini telescopes. This is part of a long term 
study to obtain both imaging and spectra for these G CS. Our first 



spectro scopic resul ts were published for NGC 3379 ( Pierce et al 



2006a[) NGC 3923 jNorris et al.l2008l) and NGC 4649 



Pierce et al 



20061: : lBridges et alj20od) . In particular, we aim at a detailed char- 



acterization of the most relevant photometric features of each GCS, 
such as colour bimodality, galactocentric colour gradients, spatial 
distributions, and the presence (or absence) of the so called "blue 
tilt" in the colour-magnitude diagrams of the GCS. All of these 
results will be later analyzed in combination with the final spec- 
troscopic sample. Additionally, a search for Ultra Compact Dwarf 
(UCD) candidates was performed in each galaxy. 

The paper is organized as follows. The galaxy sample is de- 
scribed in section 2, photometric observations and data handling in 
section 3, results in terms of colour-magnitude and colour-colour 
diagrams, colour histograms, spatial distributions, galactocentric 
colour gradients and GC integrated luminosity functions are pre- 
sented in section 4 and, finally, a summary is given in section 5. 



2 THE GALAXY SAMPLE 

Our sample consists of five nearby early-type galaxies that cover a 
range of mass. These galaxies have been imaged in several filters 
over multiple pointings with the GMOS instrument on the Gemini 
telescopes. Basic data for the sample studied in this paper, ordered 
by decreasing luminosity (and presumably decreasing mass) are 
listed in Table[T] The sample comprises two SO and three ellipticals 
(E), with B band luminosities in the range Mb = -19.94 to -21.43. 
They are located in field (NGC 31 15), group (NGC 3923, NGC 524 
and NGC 3379) and cluster (NGC 4649) environments. In this work 
we have adopted the distance m odulus from the surf ace brightness 
fluctuation method presented bv lTonrv et al.l ( I2OO1I) . This provides 
an homogeneous source whi ch can easily be conv erted to other dis- 
tance scales, such as that of I Jensen et alj ( |2003|) . In the following 
we give a short description of each galaxy with references to earlier 
work pertaining to their GCS. 

NGC 4649 Also known as M60, this Virgo Cluster member is a 
giant elliptical (E2) within the Virgo cluster and located in a sub- 
group to the East of the main cluster concentration. NGC 4647, 
a spiral some 2.5 arcmin to the NW of NGC 4649, appears to 
be a background object since there is no evidence of reddening 
produ ced by this spiral on NGC 4649 dWhite. Keel & Conselicd 
I2OOOI) and no evidence of interaction betwee n them. The GMO S 
fields presented here, and in our earlier work dPorbes et al.ll2004l) . 



were chosen to minimize the effect of possible contamination of 
the NGC 4649 GC system by NGC 4647 GCs. 

With a luminosity of Mb = -21.43, NGC 4649 is the most 
massiv e galaxy in our sample. B ased on spectra of the central 
region, iTerievich & ForbesI J2OO2I) used Lick indices to measure 
an old age of II Gyr and a metallicity of [Fe/H] = 0.3. NGC 
4649 shows strong X-r ay emission arising from a hot gas halo 
to' Sullivan et al.l l200lh . Chandra observations reveal the pres- 
ence of numerous discrete sources that are mostly identified as 
low mass binaries (LXMBs) within NGC 4649's glo bular clus- 
ters jSarazin et a l.'2003': 'Randall. Sarazin & Irwi nl2004l) . Using the 
XMM satellite, iRandall, Sarazin & Irwin (2006) present a hot gas 
density profile of s lope -1.5 ± O.I, deter mined between radii of 
10 and 300 arc sec. iFukazawa et al. I ( I2OO6I) included this galaxy in 
their study of 53 ellipticals using archival Chandra X-ray data. Fit- 
ting a P model to the X-ray surface brightness profile they obtained 
a value of /3=0.512, corresponding to a power-law slope of -2.07. 

The NGC 4649 GC S was studied by 
ICouture. Harris & AUwrightl ( Il99lh in the B and V bands. 
Within the small CCD (2.1 x 3.4 arcmin) area, no bimodality was 
detected and the authors derived a mean colour of (B-V) = 0.75 
based on 82 globular cluster candidates. 

Bimodality in the GCs colour distribution was detected by 
iLarsen et al.] ( I2OOII) on the basis of HSTAVFPC2 observations that 
found colour peaks at (V-I)o = .95 and 1.26. This bimodality was 
confirmed bv lForbes et al. 1( |2004 in the g' and i' bands using Gem- 
ini/GMOS imaging. These authors estimated a total GC population 
of 3700 ± 900, giving a specific frequency Sn= 4. 1 ± 1 . Th e se val- 



ues are in good agreement with those presente d bvlLee et al. (200S ) 
in their Washington CTi photometric study. IPierce et al.. C2006b ) 
found that most of their 38 spectroscopically-confirmed GCs are 
old (> 10 Gyr), but they found some y oung and super-solar metal- 
licity GCs. This sample was used by iBridg eTet al. (2006) to an- 
alyze the GC kinematics and dark matter content of NGC 4649 
within 3.5 effective radii (20 kpc). No rotation of the GCS was 
detected, and both isotropic and axisymmetric orbit-based models 
gave support for the presence of a dark matter halo in NGC 4649. 



NGC 3923 This is an E4-5 galaxy within a small group and ex- 
hibits a well known shell structure dMalin & Carter Il980l) . These 



kind of structures are freque nt in ellipticals and SO ga laxies not 
belonging to galaxy clusters dSchweizer & Seitzeilll992l) . and they 
are usually identi fied as the resu lt of a n interaction or merger (e.g. 
[Thomson & Wrig ht 1990; Ouin nll984l) . These shell structures may 
be long-lived and therefore an interaction/merger event may have 
occurred long time ago. 

From sp e ctrosc opic data, iThomas et al.l ( |2005|) and 
iDenicolo et al ] 1 I2OO5I) derive a central galaxy age of 3.4 and 
2.6 Gyrs respectively. Thus about 3 Gyrs ago, the galaxy appears 
to have experienced a central starburst (perhaps induced by the 
same eve nt that created th e shells ). The galaxy has a halo of hot 
gas (e.g. lO'Sullivaiietail ( l200lh ). Using the ROSAT satellite, 
i Buote & CanizaresI h99% fit a /3 model to the X-ray surface 
density profile out to 500 arcsec and d erived a power-law slope of 
-1.76 ± 0.02. iFukazawa erakl ( I2OO6I) fitted a double /? model to 
the X-ray surface brightness and found /3 values of 0.598/0.314 
for the less/more extended component. These values correspond to 
power-law slopes of -2.54/-0.88. 

Two previous photometric studi es dZepf. Geisler & Ashmanl 
ll994l : IZeDf, Ashman & Geisleilll995l) using ground-based Wash- 
ington C and Ti photometry, noted the presence of colour peaks 
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Table 1. G alaxy sam ple. Coordinates and Hubbl e types are t aken from the NASA Extragala ctic Database, B magnitu des from the RC3 catalogue 
ide Vaucoule urs et al. 1 99l| ), distance modulus from lTonrv et alj j200ll) . extinction from Schlegel, Finkbeiner & DavisI mm andX -ray luminosity is from 
lO'SuUivan et al.U2Q0ll) . 
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at (C-Tl) = 1.47 and 1.87, somewhat redder than typical values. 
This led these authors to suggest a higher than average metal con- 
tent for the NGC 3923 GCs. The total globular cluster population 
was estimated to ran ge from 2000 to 7000 and the Sn, from 4.4 
to 7.3. More recentlv lsikkema et all l l2006h using HST/ACS imag- 
ing measured similar blue and red GC projected density slopes of 
-0.87. However, their fits were obtained within the central 100 arc- 
sec and could be affected by the flattening of the inner GC profiles 
as seen in many other galaxies. They obtained a Sjv value of 5.6, 
and concluded that NGC 3923 has the highest Sjv of any isolated 
ellip tical. 

iNorris et al.l ( l2008h found that the GCs spectroscopically ex- 
aminated are old (> 10 Gys) with [Z/H] = -1.8 to 0.35. Addition- 
ally, they show that the diffuse light of the galaxy at 3Re is found to 
have ages, metallicities and [a/Fe] abundance ratios indistinguish- 
able from those of the red GCs. 

There are two spectroscopically confi rmed UCDs in this 
galaxy, found bv lNorris & Kannappanl ilOl ih . 



NGC 524 This galaxy is part of a gr oup that includes at leas t eight 
smaller galaxies known as CfA 13 jGeller & Huchralll983 ). Thi s 



group was detected in the X-ray study by iMulchaev et al. N2003h . 
The P parameter fitted to this galaxy is 0.41 ±0.02, which translates 
into a power-law slope of -1.46±0.12. The galaxy nucleus appears 
to be chemically decoupled from its bulge, being about 0.5-0.6 dex 
more me tal-rich according to Isil'chenko, Afanasiev & VlasvukI 
1 19921 and lSil'chenkol ilOOd) . who also find evidence for a dusty 
disk and ionized gas within 3 kpc of the galaxy centre. 

[Harris & HanesI ( fl98^ found a rich and extended GCS with 
some 2830 ± 880 clusters. iLarsen et al.l ( 1200 ll) used a KMM anal- 
ysis on HSTAVFPC2 dat a to find colour peaks at fV-ZJo = 0.98 and 
1.19 mag. More recentlv. iBeaslev et al. I ( |2004|) found a broad (V-I) 
colour distribution and the possible existence of peaks at (V-I)a = 
0.9 and 1.1. However, none of these results are conclusive. Beasley 
et al.'s Keck/LRIS spectra of a small sample of GCs points to a 
metallicity range from [Fe/H] = -2.0 to 0.0 and a decreasing trend 
of the [a/Fe] ratio with metallicity. The GCs are generally consis- 
tent with being old. Regarding UCDs, there is one confirmed object 
by Norris & Kannappan (private communication). 

NGC 3115 This is a relatively isolated SO galaxy with a domi- 
nant bulge, located to the South of the Leo group. The X-ray emis- 
sion from NGC 31 15 is low, co ntained within 10 kpc, and probably 
dominated by discrete sources jFukazawa et al 1 l2006h . Its GC sys- 
tem is bimodal with colour peaks at {V-1)q = 0.96 and 1.17, and has 



a low local Sjy value of '--^ 1.3 within the inner region of the galaxy 
jKundu & Whitmor3l 19981) . These authors also find that while the 
red GCs subpopulation has a spatial distribution comparable to that 
of the thick disk, the blue GCs exhibit a distribution more similar 
to the bulge/halo. 

Selecting from the iKundu & Whitmord i ll 9981) GCs sample, 
iKuntschner et a D J2003) derived spectroscopic ages and metallici- 
ties for 17 GCs, finding both subpopulations to be ~12 Gyr old and 
coeval within 2 to 3 Gyr. This result is consistent with lPuzia et al.l 
( I2OO2I) who also find coeval ages within 3 Gyr. 

From GMOS spectro sc opy of the galaxy halo, 
iNorris, Sharpies & Kuntschned (l2006h calculated an age of 
10 to 12 Gyr, similar to that of the GCs. They also obtain an [a/Fe] 
ratio of ~ 0.2-0.3, consistent with the typical GC values. 

NGC 3379 This galaxy is the domi nant object of the nearby 
Leo group. iTerlevich & ForbesI (|2002[) derive an old central age 
of ~9.3 Gyr, consistent with ages larger than 8 Gyr deter- 
mined thro u gh HS T-NICMOS observations of resolved stars by 
Gre gg et al.. (|2004|). The g a laxy c ontains a very modest hot gas 
halo, and iRevnivtsev et aU ( I2OO8I) suggest that a significant con- 
tribution to the observed X-ray flux of the galaxy comes from un- 
resolved point sources. 

From a CCD study of the NGC 3379 GC svstem. lAihar et all 

j 19941) detect ed ~ 60 GCs but no evidence of bimodality. Using 
HSTAVFPC2. lLarsen et all bOOll) confirmed that NGC 3379 con- 
tains a poor GCS. B oth IWhitlock, Forbes & BeaslevI | |2003|) and 
iRhode & Zepj l l2004h found evidence of bimodality using broad 
band colours. These last authors estimated a total GC population of 

270 clusters giving a low Sjv valu e of 1.2 ± 0.3 . 

In their spectroscopic study, IPierce et"aL ] ( l2006al) found the 
GCs to be consistent with old ages, i.e. 10 Gyr, with a wide range 
of metallicities. 



3 DATA 

3.1 Observations and Data Reduction 



The images presented here were taken using GMOS dHook et al.l 
l2004lfl in imaging mode on both the Gemini South and Gemini 



^ Further information on GMOS can be found on the instrument homepage 
jhttp://www.gemini.edu/sciops/instruments/gmos/) 
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Figure 1. R-band images of the sample galaxies from the DSS2. Each im- 
age has a field-of-view of 15 X 15 arcmin. From top to bottom we show 
NGC 4649, NGC 3923, NGC 524, NGC 3115 and NGC 3379. For each 
galaxy, we display the red GC candidates in the right panels, and the blue 
GC candidates in the left panels. The GMOS fields are also depicted. North 
is up and East, to the left. 
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North telescopes. The instruments consist of three 2048 x 4608 
pixel CCDs, with a scale of 0.0727 arcsec pixel"^. We used 2x2 
binning, yielding a 0.146 arcsec pixel" ^ scale. The GMOS cam- 
era has a field-of-view of 5.5' x 5.5'. Except for NGC 524 and 
our Comparison Field (where we had to reject and repeat some 
bad images), four images per field were taken for each of the three 
SDSS filters g', / and i' l lFukugita et 19961) . The telescope was 
dithered between exposures to facilitate cosmic-ray removal and to 
fill the gaps between the CCD chips. Details of the observational 
parameters are given in Table |2] In addition, we make use of the 
William Herschel Deep Field to help estimate contamination (see 
Sec 13.41 ). Fig. [T] shows the positions of the GMOS fields around 
each galaxy superimposed on Digital Sky Survey images. 

With the aim of obtaining a homogeneous photometric sam- 
ple, we re-reduced the images provided by Gemini. The raw images 
were processed using the Gemini GMOS package within IRAfQ 
(GPREPARE, GBIAS, GIFLAT, GIREDUCE, GMOSAIC). In all 
cases we chose suitable sets of raw bias and flat-field frames from 
the available data in the Gemini Science Archive (GSA). The im- 
ages taken with the i' band filter on Gemini South show consider- 
able fringing. We therefore built suitable fringe images from seven 
exposures of blank sky regions taken as part of the baseline calibra- 
tion the night after the observations of NGC 3923 and NGC31 15. 
These frames were used to correct /' images (listed in table O by 
means of the GIFRINGE and GIRMFRINGE tasks. The resulting 
images for each filter were then co-added and cosmic-ray cleaned 
using IMCOADD. Because the images were taken in queue mode, 
some of them were obtained under slightly different sky conditions. 
As a consequence, we set "fl.scale = yes" in order to combine 
images with different FWHM and signal. In this step, we tuned 
some IMCOADD parameters such as "aperture", "statsec", and set 
"scalenoise ~ 3" in order to avoid losing good pixels and spurious 
zero points offsets during the adding process. These final co-added 
images were then used for all the subsequent data analysis. 



3.2 Photometry 

Galaxy light subtraction was done using an iterati ve combination of 
the S Extractor background modelling approach dBertin & AmoutsI 
19961) and median filtering, in a similar way to that presented in 



Puzia et alj j2004h . This procedure was implemented in a script 



which yields a catalogue of all the objects detected by SExtrac- 
tor and a galaxy light subtracted image. This software gives bet- 
ter res ults than IRAF tasks such as DAOFIND. As in Fo rbes et al.l 
( |2004|) . the ;' band objects list was adopted as input for the sub- 
sequent photometry. At the distance of the targets listed in Table 
[T]most of the GCs would be ex pected to be u nresolved (starlike) 
sources. The Daophot package l lStetson|[l987l) within IRAF was 
used to obtain psf magnitudes for all objects detected by SExtractor, 
and they were then separated into resolved and unresolved objects 
following the outline presented in Section 1331 After that, aperture 
photometry was obtained for all the resolved objects. For this pur- 
pose we have used the PHOT task with a fixed 2 arcsecond aper- 
ture. This means that we have good photometry for unresolved and 
marginally resolved objects. A "master" catalogue of resolved and 



^ IRAF is distributed by the National Optical Astronomical Observato- 
lies, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation 
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Table 2. Summary of observations. This table lists the program ID, numbers 
of the fields, exposures and final seeing of the co-added images. 



Galaxy 


Gemini ID 


Fields 


Texp.(s) 


FWHM (arcsec) 








g' r' •' 


/ / •/ 

g r 1 



NGC 4649 


GN-2007A-Q-37 


1 


4x120 4x120 


4x120 


048 


0.49 


0.48 






2 


4x120 4x120 


4x120 


048 


0.45 


0.45 




GS-2007A-O-49 


3 


4x120 4x120 


4x120 


1.11 


1.05 


1.12 




GN-2007A-Q-37 


4 


4x120 4x120 


4x120 


0.51 


0.51 


0.47 


NGC 524 


GN-2002B-Q-25 


1 


4x150 4x100 


4x100 


0.67 


0.60 


0.65 






2 


6x150 4x100 


5x100 


0.85 


0.76 


0.72 






3 


4x100 4x100 


5x100 


0.79 


0.81 


0.84 


NGC 3379 


GN-2003A-O-22 


1 


4x200 4x100 


4x100 


1.03 


0.99 


0.84 






2 


4x200 4x100 


4x100 


0.54 


0.54 


0.52 






3 


4x200 4x100 


4x100 


0.93 


0.98. 


0.90 


NGC 3115 


GS-2004A-Q-9 


1 


4x200 4x100 


4x100 


0.77 


0.65 


0.53 






2 


4x200 4x100 


4x100 


0.80 


0.84 


0.70 


NGC 3923 


GS-2004A-Q-9 


1 


4x200 4x100 


4x100 


0.79 


0.68 


0.67 






2 


4x200 4x100 


4x100 


0.72 


0.62 


0.55 






3 


4x200 4x100 


4x100 


0.72 


0.63 


0.59 


Comp. Field 


GN-2003A-Q-22 


1 


3x200 3x100 


3x100 


0.64 


0.55 


0.50 


WHDF 


GN-2001B-SV-104 


1 


6x300 6x300 


6x300 


0.86 


0.67 


0.67 


Blank-sky 


gs-cal20040120 


1 




7x300 






0.47 



unresolved objects was built by combining photometry of objects 
successfully measured in all three filters. 

We searched in the GSA for standard star fields observed dur- 
ing the same nights as our targets and the Comparison Field. They 
were re-reduced using the same biases and flats applied to our 
scienc e frames. The synthetic transformation from iFukugita et aU 
( 1 19961) was adopted to obtain the g , r , i standard magnitudes. 
Because of the small number of the available standards, usually 
3-4 stars, and their very limited range in airmass, this calibration 
was suitable only for establishing a mean zero-point value for our 
frames, and not for a more complete evaluation of the coefficients 
in airmass or colour (see Table |3). Thus, we measured these zero 
points, and we used the atmospheric extinction coefficients given 
by the Gemini web page. 

NGC 4649 is an exception because standards were not ob- 
served at the same night as our program fields. So, in this case, we 
adopted the mean zero points obtained from the three sets of stan- 
dards (about five stars per field) observed during the Gemini North 
CMOS run of March 2007. 

For each galaxy, the field with the best overall seeing in the 
three filters and having a coiTesponding standard stars field, was 
adopted as the reference for the photometry. A list of common 
objects in each field was used to obtain any small zero point dif- 
ferences among the pointing s. Finally, we applied the galactic ex - 
tinction coefficients given bv lSchlegel. Finkbeiner & DavisI ( fl99^ , 
and used their Table 6 to transform them into Agi, A^' and A;/ . Thus 
we quote extinction corrected magnitudes and colours in this paper. 



3.3 Object classification and completeness 

The classification between resolved and unr esolved sour c es wa s 
carried out using the same procedure as in iForbes et al.l (|2004|) . 
Briefly, we used a combination of SExtractor Stellarity Index (0 
for resolved objects and I for unresolved ones), and the posi- 
tions of each object in the aperture mi nus PSF magnitude vs. PSF 
magnitude diagram jForte et al]l200lL and see Fig. 1 and 2 from 
iForbes et al.l2004l ). This a pproach has been tested using our follow- 
up spectroscopic studies jPierce et al. I l2006al; [Bridges et alj|2006l ; 
iNorris et al.ll2008b , and combined with three-colour selection, has 



been very successful in giving a low contamination rate by back- 
ground galaxies. 

In order to quantify the detection limits of our photometry a 
series of completeness tests were carried out on each field using the 
ADDSTARS task. The i' band magnitude range from 18.5 to 26.5 
mag was divided into intervals of O.I mag. Then, trying to avoid 
crowding effects, 200 point sources were added to the original i' 
images in each magnitude bin and the same procedure of galaxy 
light subtraction, detection and classification was followed to re- 
cover and classify the added sources. The positions of the artificial 
objects were created using a stochastic-generating Fortran program 
that follows the slope of the GC system spatial distribution. A total 
of 32000 artificial objects were added in two independent experi- 
ments per field. 

As an example. Fig. |2] shows a typical result of this procedure. 
In this case, the results for each of the i' fields of NGC 3923 are 
plotted. The left panel shows the fraction of artificial objects re- 
covered versus the input magnitude for each field. The right panel 
shows the completeness as a function of Rgai for three different 
magnitude bins in field 1. If a global completeness level is selected 
from the left plot, the right one shows that most of the objects are 
lost at small Kgai- That is, the completeness has a strong spatial 
dependence in the sense that the loss of objects is always larger at 
small radii. Taking this into account, and in order to ensure com- 
pleteness levels greater than 50% and 90% for Rgai > 45 arcsecs, 
we selected the /' band limits listed in Table[3] The 50% level was 
used to define a sample of GC candidates, but we used the 90% 
level in the subsequent analysis. 



3.4 Comparison Fields 

In addition to the galaxy sample listed in Table [T] we include 
two co mparison fields. One of them is the same as in lForbes et al.l 
( l2004h . However, here we present re-reduced and re-calibrated data 
for this field. 

We have applied the same object detection and classification 
procedures as for our program images. Fig. |3] shows, in the up- 
per panel, the colour-magnitude diagrams for the unresolved and 
resolved objects in this field. For um'esolved objects, we used psf 
magnitudes, while for resolved objects we used aperture photome- 
try with 2 arcsec apertures. Most of the objects with magnitudes i'o 
< 23.5 and colours within the ranges adopted for GC candidates 
are resolved and were correctly classified (see Sec 14. 11 . However, 
for i'o ~ 23.5 the situation is very different and there are many point 
sources displaying (g'-i') colours similar to those of GCs. Most of 
them are probably distant background galaxies. 

In order to test the results obtained with our Comparison Field 
(CF), we applied ou r classification recipe on the William Herschel 
Deep Field (WHDF, lMetcalfe et alj|200lh . This field was observed 
with the same instrument and filters as our CF. The seeing in both 
fields is similar, so the only difference are the total exposure times 
(see Table[2j. The 90% and 50% completeness levels obtained for 
our comparison field and the WHDF are listed in Table |3] 

Fig. [3] shows, in the bottom panel, the colour-magnitude dia- 
grams for unresolved and resolved objects detected in the WHDF. 
Clearly, there are two main differences between the two CMDs in 
Fig[3] First, the numbers of star-like objects brighter than !o< 23.5 
present in these fields are different. There are more stars in our CF 
than in the WHDF. This is easily seen in the red side of the CMD 
of the unresolved objects {(g'-i' ) > 1.5). 

The other difference is the appearance of the CMD at the low 
brightness extreme. The number of unresolved objects in the CF 
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Figure 2. Completeness factors for the i images of NGC 3923. The solid, dotted and dashed lines correspond to fields 1, 2 and 3, respectively. Left panel: 
Fraction of recovered objects versus the input magnitudes. Field 1 (centered in the galaxy) shows a smaller overall completeness than fields 2 and 3. The solid 
lines are the 90% and 50% levels. Right panel: Completeness for field 1 versus galactocentric radius for three different magnitude bins. The fraction of lost 
objects is strongly dependent on galactocentric radius for R<100 arcsec. 



Table 3. Ranges adopted for GC candidates selections, zero points, and completeness levels in the Jq band. 



Galaxy 


»0 


(9' - 






(r' - i')o 


Zp g' 


Zp r' 


Zp i' 


Completeness 




















90% 


50% 


NGC 4649 


19.5-24.3 


0.40- 


1.45 


0.35-0.95 


0.0-0.60 


27.885±0.005 


28.268±0.004 


28.190±0.006 


23.90 


24.30 


NGC 3923 


20.0-23.7 


0.40- 


1.40 


0.30-0.90 


0.0-0.60 


28.239±0.006 


28.128±0.008 


27.7 11 ±0.009 


23.30 


23.70 


NGC 524 


20.4-23.5 


0.40- 


1.40 


0.30-0.90 


0.0-0.60 


27.630±0.018 


27.992±0.028 


27.777±0.030 


23.00 


23.50 


NGC 3115 


18.4-23.5 


0.40- 


1.40 


0.30-0.90 


0.0-0.60 


28.239±0.006 


28.128±0.008 


27.711 ±0.009 


23.00 


23.50 


NGC 3379 


18.6-23.2 


0.40- 


1.40 


0.30-0.90 


0.0-0.60 


27.985±0.006 


28.215±0.005 


27.961±0.004 


22.90 


23.20 


Comp. Field 












27.905±0.017 


28.189±0.028 


27.977±0.017 


24.00 


24.33 


WHDF 












27.934±0.020 


28.184±0.029 


27.942±0.030 


24.75 


25.05 



drops steeply at iq ~ 23.5, whereas in the WHDF it does not hap- 
pen until !q ~ 24.2. The reason for this behaviour is that the pro- 
cess of classification improves with better S/N (i.e. with the deeper 
WHDF exposures we can make a better classification), and that the 
stellar content is not the same in both fields, as is to be expected 
given their different galactic latitudes. However, we note that the 
number of stars with colours adopted here for GC candidates is only 
30% larger in the CF than in the WHDF. Comparing the star counts 
in this field with the number from the Besancon Galaxy models 
( iRobin et alj2003h we found a very good agreement for iq < 23.5. 
For fainter magnitudes the Besancon star counts are always smaller 
than the observed counts, possibly indicating a veiy strong contam- 
ination by unresolved background galaxies. 

The analysis given in lFukugita et alj ( Il995h suggests that low 
redshift S, SO and E galaxies and medium redshift Irr galaxies can 
show colours in the ranges adopted for the GC candidates. How- 
ever, under good imaging conditions (i.e., sub-arcsecond FWHM) 
and making a careful image analysis, most of them are coiTectly 
classified as extended sources. The medium and high redshift 
galaxies are difficult to resolve from ground-based observation. 
However, most of these objects can be rejected by using (g'-i'), 
(g'-r) and (r'-i') colours. 

This means that, for iq< 23.5, the nature of the contamination 



is mainly stellar. At fainter magnitudes, the contamination arising 
from faint background galaxies increases steeply. 

In the subsequent analysis, we adopt an average of the counts 
in both fields to obtain the contamination levels of our photometry 
in NGC 3923, NGC 524, NGC 3 115 and NGC 3379. In the case of 
NGC 4649, we chose to use only the WHDF as a comparison field 
for two reasons. The first is that the photometry in the CF is not 
deep enough. The second reason is that, as we will see in Section|4] 
the colour-magnitude diagrams of the objects detected in our four 
NGC 4649 fields show a very small number of objects outside the 
regions occupied by GC candidates. The main difference between 
the WHDF and our CF is the red side of the colour distribution. 



4 RESULTS 

4.1 Color-magnitude and colour-colour diagrams 

The colour-magnitude (CMDs) and colour-colour diagrams are 
presented in Fig. |4] for all point sources detected in our GMOS 
fields (small dots), ordered by decreasing Bt brightness of the host 
galaxy. Magnitudes and colours were corrected for galactic extinc- 
tion as indicated in Section [T2l The mean (g'-i') photometric errors 
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the 50% completeness level in the case of our photometric sam- 
ple definition, and by the 90% completeness level for the analysis 
sample (that was used in the analysis presented in the following 
Sections). 

For the high brightness end, the cuts were chosen in order to 
include all GC candidates with absolute magnitudes fainter than Mi 
~ —12. The most massive galactic GC-like object, uj Cen, has an 
integrated magnitude of Mi ~ —11. However, MSl's most massive 
GC, Gl, is considerably more luminous. Also, this up per limit is in 
agreem ent with the Mv ~ —11 value suggested by iMieske et al.1 
( l2006al) t o separate bright GCs fr om UCDs, and is similar to that 
found bv lWehner & HarrisI ( l2007h for the Hydra UCD candidates. 
Adopting this upper limit, and util izing common objec ts between 
our GMOS photometry and that of iLarsen etaTI bOOll) from HST 
data, we obtain the (g' ,i' )-(V,I) transformations: 

7o = id -0.491(±0.003) (1) 
Vo = g'o - 0.34(±0.02)(g' - i')o - 0.03(±0.02) (2) 



I .... I .... I .... I .... I 
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Figure 3. Photometric properties of the sources detected in our comparison 
field (top panel) and in the WHDF (lower panel). In each panel, the left/right 
plot is showing the colour-magnitude diagram for unresolved/resolved 
sources. Mean colour en'ors are shown by the horizontal bars in each panel. 
It is possible to efficiently separate background galaxies from Milky Way 
(MW) stars only for objects brighter than ig ~ 23.5. 



are shown as small bars at (s'-i')o~ 0. The dotted and long-dashed 
lines indicate the 50% and 90% completeness limits from our com- 
pleteness tests, respectively. 

In all diagrams the GCs can be easily seen as a group of ob- 
jects clustered around (g'-i')o ~ 0.8, (g'-r')o ~ 0.6 and (r'-i')o 
~ 0.3. In order to obtain a clean sample of GC candidates (reject- 
ing MW stars and unresolved background galaxies), but allowing 
a reasonable range of GCs colours, we adopted the limits listed in 
Table[3] Regarding the magnitude limits, cuts were applied only on 
the ;■' band magnitudes. The low brightness ends were defined by 



(V - /)□ = 0.80(±0.02)(g' - i')o + 0.32(±0.02) (3) 

Taking into account the distance modulus from Table [T] we obtain 
the upper cut values listed in Table [3] which roughly correspond 
to Mi = —12 in each galaxy. Fig. |4] shows as filled circles all the 
point sources falling in the colour and magnitude ranges listed in 
Table[3]which we take as our GC candidates. 

Most of the unresolved objects in our sample define clear se- 
quences consistent with the locus of MW stars. The sparser appear- 
ance of the colour diagrams corresponding to NGC 524 is probably 
due to the slightly lower quality of the photometry. In the case of 
NGC 4649, we have higher quality and deeper photometry. The 
50% completeness level is fainter than that in any other galaxy and 
therefore we expect a slightly higher level of contamination by faint 
background objects (as shown in Section [J!4l the number of back- 
ground objects grows steeply for 23.5). 

In both colour-colour diagrams the GC candidates form a short 
sequence which merges with MW stars in the blue extreme. Thus, a 
small fraction of the candidates with extreme colours (g'-i')o^ 0.5, 
could in fact be Galactic stars. 

As mentioned before, at the distances of the targets listed in 
Table [T] we can expect that the normal GCs will look like unre- 
solved objects. How ever , in the case of NGC 4649 and NGC 3923, 
iBridges et all ( |2006|) and lNorris etaTI jloT ij) have found marginally 
resolved objects with colours similar to blue GCs which were con- 
firmed as members of the NGC 4649 and NGC 3923 GC systems. 
Considering that most of the images in Table [3] were taken with 
sub-arcsecond image quality, and because two of our targets are rel- 
atively nearby systems, we have checked the possibility that other 
marginally resolved GC-like objects exist in the lists of extended 
ones. Therefore, we inspected by eye all objects classified as re- 
solved, but having SExtractor shape parameters indicating round- 
ness similar to that of point sources in each image. As noted in Sec- 
tion l3.3l good S/N is necessary to robustly determine these param- 
eters, therefore we only analyzed objects brighter than ;q= 22 — 23 
mag. Any object showing a complex structure was rejected and 
aperture photometry was performed on the remaining ones. They 
are indicated as blue open circles in the CMDs, and with the same 
symbols as GC candidates in the colour-colour diagrams. This sam- 
ple of marginally resolved objects is expected to have a higher level 
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of contamination by background galaxies than our final sample of 
GC candidates. 

Since th e half-light radii o f UCD candidates (5-100 pc, 
according to iMieske et alj |2008|) are larger than those typical 
for GCs, we have extended the inspection described in the last 
paragraph to higher luminosities, including both unresolved 
sources and marginally resolved ones. Plotting the point source 
UCD candidates in the colour-colour diagrams (see (g-r) vs. 
(r'-i') panels for NGC 4649 and NGC 3923 in Fig. EJ, they 
seem to define two distinct groups. One of them includes objects 
falling exactly on the GCs sequence. The other one falls on the 
MW stellar sequence, which merges with that of the GCs at the 
blue extreme. Considering this, we have classified our unresolved 
UCD candidates as follows: "type I" if they belong to the first 
group (the most secure, shown as red filled squares in Fig. |4), 
"type III" if they belong to the second group (less secure, and 
probably MW stars, shown as green filled pentagons in Fig. |4). 
Additionally we have defined "type H" as those objects that seem 
to overlap both groups (red circles in Fig. |4]. Some marginally 
resolved UCD candidates, shown as blue filled triangles in Fig.|4] 
were also found. The full photometric tables are presented in the 
Appendix. In what follows we give a summary of each GCs system: 

NGC 4649 This galaxy shows a well-defined GC sequence in the 
colour-colour diagram. The bimodality in (g'-i) is clear and de- 
tectable even in the (g'-r ) vs. (r'-i') plane, where we can see two 
slightly superimposed clumps of GCs. There are bright blue and red 
GC candidates (Mj < —11), so the brightest objects suggest a non 
unimodal colour distribution. However, the most striking character- 
istic in NGC 4649 is that the bluest GCs seems to become redder 
with increasing luminosity, i.e., they show a "blue tilt" as will be 
discussed later. 

We have found 63 marginally resolved GC candidates (shown 
as blue circles), some of which were c onfirm ed as belonging to the 
NGC 4649 system by ISridees et aP ( l2006h (see Table |AT] in the 
Appendix). These objects have a broad color distribution but they 
seem to follow that of the GCs. 

Regarding UCD candidates, we have found 15 unresolved 
objects and one marginally resolved bright object. All of them 
have absolute magnitudes —14.5 < Mi < —12, well in the UCD 
luminosity range. We have 6 type I UCD candidates (including the 
one marginally resolved), 3 type II and 6 type III. 

NGC 3923 The GC bimodality is easily detectable, but only for 
objects with Mi > -11. In contrast with NGC 4649, the brightest 
objects show a broad and nearly unimodal distribution. This phe- 
nomenon has previously been observed in other giant galaxies, such 
as NGC 1399 l lOstrov et al.lll998l : lDirsch et al.l2003h . 

In this galaxy we find 12 marginally resolved objects, all 
of them with blue colours: (g'- i ')o < 1.0. A few of them have 
been confirmed by iNorris et al ij) as members of the NGC 
3923 system. In particular, object ID#285 is a clearly extended 
object. We have also found 32 unresolved objects with —15 < 
Mi < —12, (3 type I and 10 type II). The remaining 19 type 
III UCD candidates are very probably MW stars. One of our 
type I candidates, ID#336, was spectroscopically confirmed by 
iNorris & Kannapp aij ( l20Ilh as a real UCD (UCDl in their Table 
3). Regarding marginally resolved UCD candidates, we have 
found 7, a relatively large number. Two of them, ID#760 and 
ID #1030, were confirmed as membe rs of the NGC 3923 system 
by INorris etahlfcOllh and ID#243 bv lNorris & Kannappanl l l201lh 



(UCD2 in their Table 3). 

NGC 524 This galaxy shows a broad GC colour distribution with- 
out clear evidence of bimodality. This is similar to that found in 
some other GCSs like NGC 1427 jForte et al.ll200ll) . We will re- 
turn to this point in Section 1431 The GC candidates show a large 
spread in the (g'-r') vs. (r'-i') plane. We have found 8 marginally 
resolved objects with colours consistent with GCs. Curiously all of 
these objects are very red and only found in Field 1. 

Regarding UCD candidates, there are 8 candidates with 
luminosities Mi > —15. Two are type I and 6 type II. The 
classification, however, is tentative due to the large spread of 
these objects in the colour-colour diagrams. We have also found 
one marginally resolved UCD candidate. Norris & Kannappan 
(201 1, in preparation) have confirmed one UCD object which was 
not included in our sample because it is redder than our colour cuts. 

NGC 3115 This galaxy shows the most obvious bimodality in our 
sample except for the brightest GCs (M/< —10.5). Curiously, the 
colour magnitude diagram shows that the red GCs do not reach 
the same low luminosities as the blue ones do. In turn, the colour- 
colour diagram defines a very thin sequence with two GC clumps 
(as in the case of NGC 4649). 

We also find a noticeable number of resolved GC candidates. 
A total of 38 objects fall in this category, and some of them clump 
at Mi< —11. It is interesting to note th a t most of t hese objects 
are co nfirmed GCs by iKuntschner et al.l ( l2002h and IPuzia et al] 
( |2004|). There is on e resolved UCD candidate (presented as a GC 
in lPuzia et al.ll2004l) and 7 unresolved ones. One of them is type I, 
and the remaining are type III. 

NGC 3379 This is the poorest GCS and the least massive galaxy in 
our sample. However, bimodality is detectable even in the (g'-r') 
vs. (r'-i') plane. 

There are three marginally resolved GC candidates and three 
possible UCDs but all of them were classified as type III objects. 
One object, withM/~ —1 2 and and (g' - i')o^ .8 (which was spec- 
troscopically confirmed bv lPierce et al.l(l2006al) as a GC), lies in the 
transition region between bright GCs and UCDs. We do not find 
any resolved UCDs. Some of the other bright GC ca ndidates are 
indeed confirmed as such bv lBergond et al.l ( l2006h and lPuzia et al.l 
( |2004|) . Therefore, this low-mass system also shows a detectable 
population of very massive GCs. 



4.2 Colour Histograms 

Raw (g'-i')o colour distributions for all GC candidates selected as 
indicated in Section 14.11 are presented in the left panel of Fig. [5] 
ordered by decreasing Bt brightness of the host galaxy. These his- 
tograms were built by counting the objects in bins of 0.06 (NGC 
4649, NGC 3923, NGC 524 and NGC 3115) and 0.08 mag (NGC 
3379), depending on the total number of candidates in each sam- 
ple (dotted lines in Fig. |5j- Additionally, a smoothed colour dis- 
tribution was created by replacing each object by a Gaussian of 
a = 0.04 mags (solid lines in Fig.|5). This a value was considered 
to be representative of the mean error in (g'-i') colours for fainter 
candidates. Each histogram is background-corrected, as described 
in the following, and the distribution of the contaminant objects are 
plotted as dashed lines. 

For the case of NGC 4649, we adopted the WHDF as our con- 
trol field. In all other cases, the average of counts in the WHDF 
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Figure 4. Colour-magnitude and colour-colour diagrams for all the point sources detected in our GMOS fields (small dots). The doted and long-dashed lines 
indicate the 50% and 90% completeness levels, respectively. The mean errors in (g'~i') are plotted as small bars in the CMDs. Big dots are our GC candidates 
selected as indicated in the text. In the colour-magnitude and in the (g'-r') vs. (r'—i'j diagrams we show the UCD candidates: type I (red filled squares), type 
II (red circles) and type III (green pentagons). The marginally resolved GCs are plotted as big blue circles and the resolved UCDs ones as blue filled triangles. 



and our Comparison Field was taken as representative of the back- 
ground and subtracted. 

As noted in Section 13.41 MW stars are the most important 
source of contamination for i[) < 23.5, while uirresolved back- 
ground galaxies are the main contaminants at fainter magnitudes. 
However, none of the targets listed in Table[T]are at very low galac- 
tic latitudes and the expected contamination level is probably less 
than 10 percent (as demonstrated by our previous spectroscopic re- 
sults). 

The left panel of Fig. \5\ shows that, except for NGC 524, all 



the GCSs appear bimodal. NGC 4649 and NGC 3115 show very 
well defined blue and red peaks of (g'-i' )o ~ 0.75 and (g'-i')o ~ 
1.0 — 1.1. NGC 3923 appears bimodal, but the separation between 
peaks is less clear. NGC 3379 shows a marked blue peak around 
(g'-i')o ~ 0.7 — 0.75, and a less conspicuous red population. 

As noted before, NGC 524 is the only target where no clear bi- 
modality is seen. However, its colour distribution looks rather broad 
and further high quality photometry will be required to clarify the 
situation. 

As shown in Fig. [T] the areal coverages are different in each 
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Figure 4. Continued. 



observed GCS. As a consequence of this, we have attempted to re- 
cover the global colour distribution of the entire GCS by correcting 
for areal incompleteness. In order to do this, we split each GCS into 
different galactocentric radial bins. As seen in Section |373l the inner 
regions always have lower photometric completeness. We therefore 
exclude the zones within Rgai < 45 arcsecs. The outer zones have 
very small areal completeness (and hence introduce noise), so the 
points with Rgaj > Vmax were excluded too. Here rmax was taken 
as the Rga; which shows areal completeness smaller than 20% (in 
Table |4] we listed the values of Tmax for each GCS). After that, 
we calculated the fraction of each ring "effectively" observed. And 
then, the inverse of these fractions were used as a correction fac- 
tor The number of radial bins was three in the two most massive 
galaxies (NGC 3923 and NGC 4649), and two in the other less 
populated ones (NGC 3115, NGC 3379 and NGC 524). Finally, the 
total colour distribution was obtained by adding the colour distri- 
butions of each individual radial bin. 

In the right panel of Fig.|5]we show the resulting background 
corrected histograms and smoothed colours distributions. The areal 
corrections lead to more prominent blue GCs populations. This re- 
sults from the fact that blue GCs are less concentrated towards the 
center of the host galaxy than the red ones. 

In the case of NGC 3923, the colour-magnitude diagram 
shows a unimodal colour distribution for the brightest GCs. Then 
we have split the sample in GCs having M/< —11 (long dashed 
line in Fig.|5j and Mi> —11 (the same codes as the other GCSs). 
Doing that, we see that the bright GC candidates are not bimodal. 



For NGC 3115 and NGC 3379, the small number of bright GCs 
prevents a similar analysis. 

In order to characterize the presence of bimodality in the 
integrated GC col ours, we have run RMIX on our histograms 
l lMacDonalj|2007l) . This program allows the user to choose the 
kind of distribution, number of components, and constraints. In this 
work we adopted two normal components, allowing RMIX to fit the 
position of the peak, the standard deviation and the number of GC 
candidates in each populations. Using reasonable starting values, in 
all cases RMIX converged easily, so we did not constrain any pa- 
rameter. The results are shown in Table |4] As a comparison we list 
the modal values for all the GCSs showing two clear peaks in their 
colour distribution. In the case of NGC 4649, we list two values 
for each parameter The first set of values are obtained using the 
whole sample, and the second one taking GC candidates with 21.5 
< i'o < 23.9. In this last sample both peaks are bluer In particular 
the blue peak moves 0.025 mag as a consequence of the "blue tilt" 
phenomenon. 

In the case of NGC 524 we specifically tested different kinds 
of fits. We found that a double Gaussian distributions fit was 
strongly preferred over only one leading to a x value 2.7 times 
smaller This supports the conclusion that NGC 524 is also a bi- 
modal system. However, the blue peak in this galaxy {(g'-i')o = 
0.815 ± 0.021) is the reddest one in our sample. Deeper photome- 
try could help to obtain a more conclusive result for this galaxy. 

Fig.|6]shows the positions of the GC colour peaks for the dif- 
ferent cluster populations as listed in Table|4]versus the galaxy Mb . 
Blue and red peaks are depicted as filled triangles and squares, re- 
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Figure 5. Left panel: Dotted lines are the (g'-i' )o colour histograms of all 
our GC candidates selected from samples with 90% completeness level. 
The solid lines are smoothed colour distributions, obtained as described 
in the text. The backgrounds used to connect the counts are shown as the 
dashed lines. Right panel: The same as left panel, but corrected for areal 
completeness as described in the text. In the case of N3923, GCs brighter 
than M/=-l 1 are show as long dashed hne in this panel. 
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Figure 6. (g'-i' )o colour for blue (filled triangles) and red peaks (filled 
squares) obtained by fitting two Gaussian to the colours histograms. We 
have included here the NG C 33 11 (open and fil led circles) and M87 (open 
and filled stars) points from lWehner et al.] bO OS!) and Harris (2009), respec- 
tively. Solid lines are the fits obtained as indicated in the text. Dashed lines 
are the results excluding NGC 3923 {Mb= -21.18) in the fits. 



spectively. One of the most striking characteristics of this figure 
is that the red peak of NGC 3923 seems to be too blue accord- 
ing to its absolute magnitude. To complement this, we have in- 
cluded the values correspo nding to NGC 3311 as open and filled 
circles jWehner et al.ll2008l) . The following (error weighted) linear 
fits were obtained for the red (eq.|4) and blue GCs peaks (eq.|4): 



i')o = -0,03(±0.01)Mfl + 0.42(±0.33) 



(4) 



(9 



-0.01(±0.01)Mi3 -H0.56(±0.22) (5) 



These fits are shown as solid lines in Fi g. [6] where we also 
include the peaks corresponding to M87 from iHairisI (l2009ah for 
comparison (open and filled stars). The M87 points are very consis- 
tent wit h our relations. However, it is important to note that lHarrisI 
l l2009ij) calibrated his M87 photometry with SDSS stars. As a re- 
sult, some mismatch m ay be expected. 

IPeng et al] l l2006h found that the slopes of {g — z) vs. Ms re- 
lations were, -0.036 and -0.026 for red and blue peaks, resp ectively. 
Using some common objects between fjordan et"al] ( |2009|) . and our 
NGC 4649 photometry, we have found that d{g — z) /d{g' - i') = 
1.36 ± 0.05. With this value, the slopes from Peng et al., translate 
to -0.026 and -0.019. Comparing with Equations ^ and ([5]l, we 
see that the red values are in very good agreement. However, our 
results for the blue peaks are consiste nt with no correlation. 

We note that IPeng et all | |2006|) obtain [Fe/H] oc l0-26±0 03 
for the red GCs which is identical to our value after adopting the 
colour-metallicity relation given in the next section. 
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Table 4. RMIX fits to the (g'-i' )q color distributions. The mean colours of the two modes /ij, (blue) and fi^ (red), their dispersions cj, and ct^, and the 
fraction of object assigned by RMIX to the blue GC population are listed. The last two columns give the galactocentiic radii of the annuli which show areal 
completeness smaller than 20% and the fraction of area effectively observed within those radii. 





Modal values 






RMIX 








(r 


Galaxy 


Mb 


Mr 


Mb 


o"b 


Mr 




h 


arcsec 




NGC 4649 


0.765 


1.095 


0.783 ± 0.006 


0.098 ± 0.004 


1.080 ± 0.008 


0.110 ± 0.004 


0.518 ± 0.024 


490 


0.41 








0.758 ± 0.007 


0.096 ± 0.004 


1.067 ± 0.008 


0.127 ±0.004 


0.442 ± 0.025 






NGC 3923 


0.758 


1.011 


0.746 ±0.011 


0.072 ± 0.008 


1.005 ± 0.023 


0.135 ± 0.014 


0.365 ± 0.069 


350 


0.58 


NGC 524 






0.815 ± 0.021 


0.086 ± 0.011 


1.033 ± 0.040 


0.076 ± 0.021 


0.713 ± 0.120 


315 


0.70 


NGC 31 15 


0.764 


1.053 


0.765 ± 0.007 


0.067 ± 0.005 


1.044 ± 0.011 


0.092 ± 0.008 


0.521 ± 0.037 


320 


0.58 


NGC 3379 


0.769 




0.765 ± 0.009 


0.036 ± 0.008 


1.010 ± 0.029 


0.135 ± 0.018 


0.408 ± 0.082 


290 


0.72 



4.3 Blue tilt 

As noted in Section [43] NGC 4649 is the only galaxy in the sam- 
ple where the blue GCs show a clear colou r-magnitude t r end, a 
feat ure called the "blu e tilt", and noticed by IStrader et al. 1 |2006i) 
and lHarris etall ilOOdt) . In contrast, all the other galaxies discussed 
in this paper show two broad color sequence s, without an obvious 
tilt (f or a further study of the NGC 3923, see lNorris & Kannappanl 
I20T l|). The GC colour distributions in these galaxies become uni- 
modal at bright magnitudes, a situation different from that in NGC 
4649, where blue and red clusters are clearly separated even at 

higher luminosities. 

In particular for NGC 4649, IStrader et all (|2006|) find a colour 
trend d((; — z)/ci2 = — 0.040±0.005 identic al, within the errors, to 
d{g - z) /dz = -0.037 ± 0.004 derived bv lMieske et"an l l2006bl) 
by combining the colour magnitude diagrams of M49, M60 and 
M87. 

In order to characterize the "blue tilt" in NGC 4649, using 
our GMOS photometry, we followed an approach similar to that 
described bv lWehneretal..(20Q8.) . GCs were grouped in magnitude 
bins containing 150 clusters each and, using RMIX, we obtained 
the positions of the blue and read peaks. This procedure delivered 
nine pairs of values depicted as blue and red dots in Fig.|7] In all the 
cases previously studied, the "blue tilt" appears to have a similar 
behaviour. Fig. |7] shows that the blue sequence seems to become 
redder at a magnitude around Mi ~ — 9 mag. In turn, and usin: 
a combined sample for six giant elliptical galaxies, iHarrisI ( 120091 
found evidence that the "blue tilt" starts to be detectable at absolute 
magnitudes brighter than Mj — —9.5 mag and exhibits a nonlinear 
form. However, in our NGC 4649 sample, the three faintest points 
of the blue and red peaks (see Fig. |7) were not well constrained 
by RMIX as they have significant colour errors. So, these points 
were not included and only those with i'< 23 were used to perform 
linear fits to the blue and red GCs groups: 

{g' - i')o = -0.032(+/ - 0.003)io + 1.478(+/ - 0.072) (6) 



{9' - i')o = -0.007(+/ - 0.006)id + 1.227(+/ - 0.1356) (7) 

These fits indicate that the r ed GCs do no t s how a significant 
trend, as noted for other galaxies ( lHarrisl2009allbt|Peng et alj2008l : 
iHarris et al.ll2010l) . 

W e stress that a second order fit, like that performed bv lHarrisI 
( l2009bl) , on our data, produces very unstable solutions. However, 
we cannot rule out that the inclusion of a better sample of faint 
objects may require a higher order fit. 



As a comparison Fig. [7] also includes the position of the blue 
and red peaks obtained by iHarrisI ( l2009ah for the M87 GCS also 
using g and i photometiy. Both galaxies exhibit remarkably similar 
blue tilts. 

The fact that the tilt is detectable in some but not all galaxies 
is an issue that s till deserves clarifi cation. A tentative inteipretation 
was presented in lForte et al.l ( l2007h who showed that a large spread 
of the characteristic metallicity scale of the blue GCs smears the 
colour-magnitude relation in such a way that the tilt becomes even- 
tually undetectable. In the specific case of the galaxies in our sam- 
ple, the situation is not so clear. Although the colour magnitude di- 
agrams for NGC 3379, NGC 31 15 and NGC 524 in Fig.|4]show the 
presence of very massive GCs (those with Mj < — 10 mag) their 
numbers are quite small thereby preventing a meaningful detection 
of the blue tilt. 

It is worth mentioning that iHarris et al. I (l2010h found a well 
defined tilt in the blue GC sequence in M 104, using a sample with 
several hundred GCs. In contrast, NGC 3923 with a comparable 
number of clusters (~ 670), exhibits no detectable tilt, suggesting 
that the size of the statistical sample cannot be the main reason 
for the absence of the tilt. We note, however, that our GCs sam- 
ple does not include a significant number of clusters brighter than 
All ~ —9.5, which conspires against the eventual detection of the 
tilt. Rather, our colour magnitude diagram suggests a "unimodal" 
distribution for these bright clusters (but see Norris and Kannap- 
pan 2011). Another question that deserves a more complete wide 
field analysis is why NGC 3923 seems to have a relatively sparse 
population of massive blue GCs. 

We note that the "blue tilt" has been identified as a possible 
instr umental artifa ct specifically affecting the HST-ACS photom- 
etry ( lKundull2008h and even dismiss ed on the basis of (V-I) pho- 
tometry of M87 jwaters et al.ll2009l) . However, the observations 
presented in this paper together with other ground based w orks 
(see also lForte et"aLll2007llWehner et al.ll20q8l lHarrisll20()9al). an d 
new HST/ACS based analyses jPeng et alj|2008l : IHarrisI l2009bh . 
strongly argue in favour of the existence of the blue tilt. 

Previous works have presented an interpretation of the "blue 
tilt" as the result of increasing met allicity with GC mass, i.e., 
a mass-metallicity rela tion (MMR). iBaihn & Harris! l l2009h and 
IStrader & Smithl(l2008l) give models where the physical interpreta- 
tion of the MMR is based on self-enrichment during a cluster's for- 
mation stage. According to these models, GCs with masses above 
10^ Mq will hold onto enough of their SN ejecta to enrich the gas 
from which stars are still being formed. A meaningful comparison 
of the results from different GCS is somewhat precluded by the 
adoption of different colour-metallicity relations. The photometry 
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presente d in this paper, combined with our previ ous s pectroscopic 
result s JPierce et all l2006al IPierce et al.] l2006bl and iNorris et all 
allows the determination of a new relation that includes 
67 globulars in three different galaxies (NGC 3923, 3379, 4649). 
Chemical abundances [Z/H] obtained vi a Lick indices , are on the 
scale defined by the synthetic models bv lThomas et alj l l2004h and 
clusters fainter than g'(,=23 were rejected due to their relatively 
large associated errors. 

The [Z/H] vs. (g'-i')o values for the calibrating GCs are plot- 
ted in Fig.[8]where a linear fit: 

[Z/H] = 3.51(±0.21)(3' - i')o - 3.91(±0.20) (8) 

is also shown. The overall trend is in agreement with a rediscussion 
of the (C-Tl )o vs [Fe/H] relation presented in Iporte et all ilOOH) 
and with oth er works where more c omplex funct ions were fitted 
(for e xample, IPeng et ^1120061 lEeee t al. 2008 and lBlakeslee et al.l 
I2OIOI) . However, it is worth emphasizing that all these last calibra- 
tions include higher order colour terms, mostly set by low metal- 
licity MW globulars, whose counterparts are scarce in Fig. |8] A 
more direct compa rison is possib l e usin g the recently published re- 
lation obtained by ISinnott et al. H2O10h . Those authors have used 
integrated colours in the USNO photometric system dSmith et al.l 
l2002h and spectroscopically measured abundances of GCs from 
NGC 5128 and the MW. Their linear (g'-i')o vs [Fe/H] relation 
can be easily translated to the f orm of Equat i on (|8) using the 
[Fe/H]=[Z/H]-0.131 relation from lMendel et all ( |200^ , leading to 
[Z/H]= 3.65 X (g' — i')o — 4.18, which compares very well with 
our results. 

The adoption of a constant (M/L) ratio for the blue clusters, 
and the above calibration, leads to 



in agreement with lMieske et al. I l2006bl who, in the particular case 
of NGC 4649, derive an exponent of 0.36 ± 0.11. 

Our colour-abundan ce relation was also app lied to the g' and /' 
photometry presented by ICockcroft et all ( 120091) for NGC 5 1 93 , IC 
4329 and NGC 331 1, leading to exponents of 0.23 ± 0.14, 0.37 ± 
0.05, and 0.47±0.03 respectively. These values are, in the average, 
30 percent smaller than the exponents obtained by Cockcroft et al. 
These differences are typical of the uncertainties connected with 
the adoption of different metallicity scales. 

4.4 Spatial distributions 

Fig. U] shows some already known systematics regarding the spatial 
distribution of both GC families, i.e., blue clusters usually exhibit 
relatively shallow distributions in contrast with those of the red 
ones, which are frequently more concentrated toward the galaxy 
centers. 

In this section we aim to quantify the spatial distributions 
through an analysis of the projected cluster density as a function 
of galactocentric distance using both a de Vaucouleurs law (r^''^) 
and a power law. We also present Sersic law fits in the case of the 
two brightest galaxies for which we have a good areal coverage. 

Blue and red GCs were separated according to a common 
"colour valley" at (g'-i'jo ~ 0.95 and the limiting ig magnitude 
was set at the 90 percent completeness level. Due to their relatively 
low flattenings, we used circular annuli with A log r = 0. 1 to 0.15 
for all galaxies except NGC 3115, for which we adopted an ellip- 
ticity of 0.55. 
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Figure 7. The "blue tilt" in NGC 4649. Blue and red points are the positions 
of the peaks of the blue and red GC subpopulations obtained as indicated in 
the text. The green circles are those for the M87 GCS from.Harris (20093). 
The red line is a linear fit to the blue points for NGC 4649. Here we see 
that the blue tilt in this galaxy appears at /q ~ 23 mag, which roughly cor- 
responds to A/j ~ -8.5 mag. 




(g'-i')o 



Figure 8. GC colour metallicity relation. The points ar e GCs from the pho- 
tometry presented in this work and metallicities from IPierce et alj|20063 . 
■Pierce et al. ,.2006bi and iNorris et al.ii2008t The solid line is a linear fit (see 
the text). 
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Figure 9. Projected radial profiles for all (filled circles), red (triangles) and blue GC candidates (squares) fit with a de Vaucouleurs-law (left panel) and with 
a power-law (right panel). The profiles were arbitrarily shifted in order to avoid overlapping. Since of NGC 3379 shows a very low number of red GCs and 
NGC 524 doesn't show a clear "valley", their GCs samples were not split into red and blue subpopulations. In the case of NGC 3115 the projected density 
profile was obtained in elliptical rings. Dashed lines indicate the radial ranges included in the fits. Open circles in the NGC 4649 and NGC 3923 figures are 
from ACS photometry (see text). 
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Figure 9. 

The resulting areal density distributions are displayed in Fig. 
[21 and the fit parameters listed in Table|5] As a general comment, we 
point out that both the de Vaucouleurs and the power law give very 
similar results in terms of the residual errors, although the former 
yields marginally better fits in the innermost regions. 

In the cases of NGC 4649 and NGC 3923, our observations 
were combined with the HST-ACS data allowing a fit to the inner- 
most regions. The GC candidates from I Jordan et al.l (l2009f) were 
used for NGC 4649, and in NGC 3923 the clusters were selected 
from our own analys is of the same ACS archive images used by 
ISikkema et al] | |2006|) . The ACS and GMOS profiles were matched 
by normalizing the counts in common selected radial annuli. 

For these galaxies, Sersic law fits in the form log cr(r) = si + 
S2 are also shown in Fig.|9] The respective si, S2 and values 
are: -0.57 ±0.19, -0.27 ±0.09, 0.36 ± 0.04 for NGC 4649, and 
-0.60 ± 0.20, -0.51 ± 0.20, 0.26 ± 0.09 for NGC 3923. 

A core-like distribution is detected in the case of NGC 4649, 
simila r to t hat found, for exam ple, in M87 by iLauer & Kormendvl 
( Il986h and lKundu et"ai] jl999l) . In turn, the Sersic law seems com- 
patible with a De Vacoucoleurs fit along the whole galactocentric 
range in NGC 3923. 
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Continued. 

In what follows we give some brief comments about the char- 
acteristics of each galaxy. 

NG C 4649. Our densit y profile improves that already pre- 
sented in iForbes et all ( |2004 due to a better determination of the 
level of the subtracted background. The slope values given in Ta- 
ble[5] in turn, a re in excellent agreement with those presented by 
lLeeetal1 ( l2008h . 

NGC 3923. Both Fig.Qand Fig.|9]show that the blue globulars 
have a very shallow distribution. In fact, the slope of the density 

Jirofile for these clusters is the lowest in our sample. ISikkema et al.l 
2006) also found very low values, although their photometry only 
includes the central region of the galaxy. 

This feature deserves some attention as it may be con- 
nected with the merger history of this galaxy as suggested by 
the well known ripp les detectable in its brightness profile (see, 
lMalin&Carteilll980h . 

NGC 524. As mentioned previously, the colour histogram of 
the GCs in this galaxy do not show a clearly detectable bimodality 
and, for this reason, the density profile displayed in Fig. |9] corre- 
sponds to t he whole cluster popul ation. The derived value is con- 
sistent with lHarris & HanesI Jl985l) . 
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We note, however, that an inspection of the innermost region 
in this galaxy, indicates that the red GCs seem v ery concentrated ( a 
situation that resembles the case of NGC 1427. lForte et al.lbOOlh . 
However, splitting the sample at {g' — i')o = 0.9 as shown by Fig. 
[T] the "reddest" GCs exhibit a clump inside Rgai < 1 arcmin. In 
contrast, this feature is absent in the "blue" clusters. 

NGC 3115. This galaxy exhibits a high flattening and its GCs 
seem to follow the brightness isophotal contours. However, our un- 
even spatial coverage does not allow a quantitative confirmation 
(e.g., via azimuthal counts). 

NGC 3379. Since this galaxy shows a very low number of GC 
candidates, the density profile displayed in Fig.|9]coiTesponds to the 
w hole GC populat i on. Th e resulting slope is steeper than that given 
bv lRhode & Zepj ( |2004|) who found -1.4 in their wide field study. 
This may indicate a difference in the background level adopted in 
these two studies. 

Fig. Q] shows that the blue and red candidates in this galaxy 
are not homogeneously distributed around the center. The red GCs 
show a clump to the NE, while the blue ones are clearly less numer- 
ous in that region. This inhomogen eous distribution is also seen in 
IWhitlock. Forbes & Beaslevl ( l2003h (their Fig. 3). 



4.5 Galactocentric colour trends 

Fig.[TO]shows the (g'-i'jo colours of all the globular cluster candi- 
dates brighter than the 90 percent completeness level as a function 
of the normalized galactocentric radius Rgai/Reff, where Reff is 
the effective radius in the B band taken from RC3 data in NED. 

Mean colours for each GC family were obtained within galac- 
tocentric bins of variable size (containing the same number of clus- 
ters) and fit with a logarithmic law 



)o ~ a log{R/R, 



eff) 



(10) 



The a and b parameters are listed in Table|6] along with the as- 
sociated uncertainties obtained through a least square fit, for each 
GC subpopulation except for NGC 3379 where, due to the low 
number of clusters, blue and red GCs were grouped together. The 
listed colour gradients were combined with Equation |8] to obtain 
values of l^[Z / H\/log{R/ R^f f) and were included in Table|6] 

The overall appearance of the Fig.[TO]suggests very mild gra- 
dients over a galactocentric ran ge larger than ~ 1 to 2 Re//. A 
similar trend has been noted by iHarrisI (l2009allbl) in M87 and six 
other GCS belonging to giant elliptical galaxies. Beyond this, each 
galaxy has its own and distinct behavior. In particular, NGC 524 
and NGC 3115 (both SO galaxies ) exhibit the highe st colour gradi- 
ents in their innermost regions ( Tortora et al .120 lOl) . 

Fig. [To] also includes the (g'-i'jo halo colours derived from 
our images. They were obtained within the central 2 arcmin of the 
galaxies because this zone is not strongly affected by the sky level 
adopted. These colours are remarkably similar to the mean colours 
of the red GCs, a similarity alrea dy noted in the inner regions of 
several galaxies (see, for example. iForbes & Fortel200ll) . 

Our data adds to the increasing evidence in the literature for 
colour, and hence metallicity, gradients in the indiv idual GC sub- 
populations (e.g. iHairisI l2009bt iForbes et al.ll201lh . Such gradi- 
ents are indicative of a dissipative formation process for both GC 
subpopulations. In a detailed study of the NGC 1407 GC system, 
IForbes et al. I1 I20T ll) suggested that a transition occurred between 
the inner region with strong GC gradients and an outer region with 
no colour gradient. This was interpreted as support for the two 



■■ ■ I - I ■ I - I' J ! ■ I ' h i' 




0.4 

1.2 



0.4 



0.4 



1.8 - 



0.4 



f ' I ' r ' I 



I ■ I ■ I ■ I 



■ M l ■ l - i I ■ I ■ I • ■ I ■ I ■ I ■ 



■ I ■ I ■ I ■ I ■ I ■ I ■ I ■ I ■ I ■ 



■. V: , '. 
0.8- ■.■.'.•.»:■. i-.v*-.'.i-....-^/-.-. 



■ I ■ I ■ I ■ I ■ I ■ I ■ I ■ I ■ 



I ■ I ■ I ■ I ■ I ■ I ■ I ■ I ■ I 



0123466789 10 

Figure 10. (g'-i' )o colour index versus projected galactocentric distance 
normalized by R^ffi in the B band). The colour profile for each galaxy 
halo, within a radius of 2 arcmin, is shown as the solid line in each panel. 
From top to bottom: NGC 4649, NGC 3923, NGC 524, NGC 3115 and 
NGC 3379. 



phase formation model of early-type galaxies dOser et al.ll2010h in 
which such galaxies are built up from a dissipative core that experi- 
ences later accretion. Our data show shallower gradients and fewer 
GCs than detected in NGC 1407 by Forbes et al. Further data on a 
rich GC system like NGC 4649 would be needed to search for such 
a transition. 

The GC colour gradients are also compared to the colour gra- 
dient of the underlying starlight in the galaxy. In each case where 
the GC system can be clearly separated into two subpopulations 
we find that the galaxy gradient matches the red GC subpopu- 
lation gradient in both absolute colour and in slope over a com- 
mon radial range. This provides strong support for the idea that the 
red GCs are associated with the bulge/spheroid stars of early-type 
alaxies, and hence may have shared a common formatio n epoch 
Forbes &Fortell200 it iForte et al.ll2009l: ISpitler et al.ll20I0l) . 
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Figure 11. Raw counts (solid histograms), and completeness corrected counts (dashed line) for the five galaxies in the sample as a function of ig magnitudes. 
The nonnalized backgrounds adopted in each case are shown as dotted lines. 
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Table 5. Slopes and errors of the surface density profiles of all, red and blue GC candidates obtained using a weighted least squares method. We list the results 
for power and de Vaucouleurs laws. 



Galaxy 




all 




red 




blue 




power law 


de Vaueouleurs 


power law 


de Vaucouleurs 


power law 


de Vaueouleurs 


NGC 4649 


-1.37 ± 0.04 


-0.67 ± 0.02 


-1.52 ± 0.04 


-0.73 ± 0.02 


-1.21 ± 0.08 


-0.58 ± 0.03 


NGC 3923 


-1.25 ± 0.07 


-0.64 ± 0.03 


-1.45 ± 0.09 


-0.73 ± 0.04 


-0.86 ± 0.08 


-0.44 ± 0.05 


NGC 524 


-1.67 ± 0.06 


-0.89 ± 0.05 










NGC 3115 


-1.59 ± 0.13 


-0.78 ± 0.07 


-1.74 ± 0.17 


-1.08 ± 0.10 


-1.45 ± 0.15 


- 0.73 ± 0.07 


NGC 3115(cllip) 


-1.97 ± 0.12 


-1.02 ±0.05 


-2.45 ±0.11 


-1.28 ± 0.08 


-1.62 ± 0.05 


-0.81 ± 0.04 


NGC 3379 


-1.81 ± 0.16 


-0.98 ± 0.08 











Table 6. Slopes and errors of the mean (g'-i')o colours in function of Rgai- These values were obtained by fitting the relation: (g' — i')o = a X 
log{Rgai/ Rftf f) + b. The [Z/H] gradients, obtained using Equation [s] are listed in columns four and live. In the cases of NGC 4649, NGC 3923 and 
NGC 3115, two values are listed. The first one show the slopes obtained by using the entire sample of points showed in the figure, and the second one, 
corresponds to a fit excluding the inner point. 



Galaxy 


slopes 


-i'la 


slopes IZ/HI 






red 


blue 


red 


blue 




NGC 4649 


-0.0I5±0.007 


-0.001 ±0.011 


-0.052 ±0.024 


-0.004 ± 0.038 


50 




-0.021 ±0.009 


0.004±0.015 


-0.073 ±0.032 


0.014 ± 0.052 


50 


NGC 3923 


-0.050±0.024 


-0.051±0.019 


-0.175 ±0.084 


-0.179 ± 0.067 


30 




-0.095 ±0.024 


-0.033±0.025 


-0.333 ±0.086 


-0.115 ± 0.088 


30 


NGC 524 


-0.099±0.043 /-0.194 ± 0.069 


-0.347 ±0.152/ 


-0.681 ± 0.245 


30 


NGC 3115 


-0.030± 0.028 


-0.076±0.016 


-0.105±0.098 


-0.266 ±0.058 


20 




-0.022± 0.037 


-0.081 ±0.024 


-0.077±0.130 


-0.284 ±0.085 


20 


NGC 3379 




-0.066 ± 0.055 




-0.232 ±0.193 


15 



4.6 GC Integrated Luminosity Functions and GC Specific 
Frequencies S jv 

4.6.1 Integrated luminosity functions 

In this section we aim to determine the integrated GC luminosity 
function (GCLF) for each galaxy. Fig.[TT]shows the raw and com- 
pleteness corrected counts of cluster candidates as well as those of 
the adopted comparison field as a function of iq magnitude. NGC 
4649 and NGC 3115 exhibit a good contrast at the expected TOM 
magnitude between clusters and contaminant objects allowing this 
kind of analysis. In NGC 3923, we don't have such a good con- 
trast. However, the completeness and the background counts are 
good enough to try to recover some useful information. In NGC 
3379 the low number of candidates, and the strong dependence of 
the background adopted make it difficult to obtain reliable results 
and the parameters of the GCLF are not very well constrained. For 
that reason the points with /' > 23.2 were rejected from the fit. In 
the case of NGC 524, we expect to have the TOM at iq ~ 23.9, 
where our photometry has a very low photometric completeness 
and therefore this galaxy was not included in what follows. 

The resulting GCLF, determined on the basis of the net counts 
within 0.25 mag bins, are displayed in Fig. [12] together with the 
normalized counts of the subtracted background. 

Both a Gaussian and a t5 function fit were performed on each 
histogram in order to determine the position of the turn-over mag- 
nitude (TOM) and the dispersion parameter As the results do not 
show clear differences between both approaches, we only give the 
parameters of the Gaussian fit in Table |7] 

Reliable estimates of the TOMs in NGC 4649 and NGC 3115 



are derived since our adopted limiting magnitude is fainter than the 
TOMs by approximately one magnitude. 

For these galaxies we also obtained the TOMs of the blue and 
red GC populations defined in terms of the (g'-i')o "valley" in the 
colour statistics. The corresponding corrected luminosity functions 
for each GCs subpopulation are shown in Fig. [13] A comparison 
shows that, in NGC 4649, the blue and red TOMs are coincident 
within the errors. This agreement is expected as the ip magnitudes 
are n ot strongly affected by metallicity effects (see lAshman et al] 
I199A . The situation is not clear in NGC 3115 where the blue TOM 
seems fainter than the red one, contrary to the expectations for two 
coeval GC populations with different metallicities. An inspection 
of Fig.[T]shows that our areal coverage of the central region of this 
galaxy is rather poor, resulting in the loss of a fraction of the red 
GCs population (because red GCs are more concentrated towards 
the center of the galaxy than the blue one) and possibly producing 
an spurious effect on the magnitude statistics of these clusters. 

For NGC 3923, and as the TOM seems very close to the mag- 
nitude cut off (at a 50 percent completeness level), we only derive 
the dispersion. In this case we adopted the I band TOM given by 
Kundu & Whitmore ( 2001) which, at the adopted distance modulus 
becomes ;o=23.8 ± 0.28. 

The adoption of the absolute TOM magnitudes given by 
iKundu & Whitmord JioOll) . and using the relation given in Equa- 
tion[T] lead to distance moduh (m-M)i= 30.92 ±0.06 and (m-M)i = 
29.75 ± 0.09 for NGC 4649, NGC 3115 respectively. 

A co mparison withlLarsen et al. I ( I2OOII) , iKuirdu & Whitmord 

( Il998l) and lKundu & Whitmord ( 1200 ll) (previous correction by the 
TOM for the MW GCs and the Mv adopted by these authors) 
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shows an excellent agreement between their distance moduli and 
ours. 



4.6.2 Specific frequency and toted GC population 

Adopting Gaussian integrated luminosity functions defined by the 
parameters given in Table|7]in combination with the projected areal 
density profiles previously discussed, we derive the cumulative spe- 
cific globular frequency Sn ~ A^totlO" *'*^^^^^' for each galaxy 
within maximum galactocentric radii of 50 and 100 kpc. The first 
radius seems a reasonable value in order to include most (if not all) 
the CGs associated with a given low mass galaxy. The largest ra- 
dius seem more appropriate fo r massive galaxies that have very e x- 
tended GCS (see, for example. iBassino et alj2006llHMi-isl2009ah . 

In order to estimate Ntot for NGC 4649 and NGC 3923, we 
have integrated the de Vaucouleurs law profiles obtained by fitting 
the GMOS+ACS profile presented in Fig. |9] In the particular case 



of NGC 31 15, we have integrated our de Vaucouleurs GMOS pro- 
file an d we adopted the WFPC2 photometry given bv lLarsen et al.! 
I2OOII to complete the counts in the innermost regions of the galaxy 
not covered by our frames. As those density profiles were obtained 
considering samples with 90 percent completeness magnitude cuts, 
these three values were corrected by calculating the fraction of the 
GCLF not included beyond these magnitude limits. 

Absolute galaxy luminosities were obtained using the photo- 
metric values given in the NED database and the distance moduli 
listed in Table[T] The total number of GCs, galaxy visual luminos- 
ity and inferred Sjv values are included in Table|7] As an indicative 
of the uncertainty in our estimation of A'^tot and Sn, we list as er- 
ror the difference between the results obtained by integrating a de 
Vaucouleurs and power law profiles. 

For NGC 4649, there are two previously published values of 
Sjv in the lite rature. Both of th em were calculated for 50 Kpc and 
they are from iLee etal] (l2008h . Sn = 3.8 ± 0.4, and that from 
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Figure 13. Corrected luminosity functions for each GCs subpopulation in galaxies NGC 4649 (left) and NGC 3115 (riglit). The LF for red GCs are shown as 
solid lines and those for blue GCs as dashed lines. 



iForbes et al. I l200i 5'jv = 4.1 ± 1. A comparison wilh Ihe num- 
bers listed in Table [7] shows that our new Sn value is very good 
agreement with those values. 

For NGC 3923 ISikkema et all |2006| have obtained Sn = 
5.6 ±1.3 and Sn ~ 8.3 ± 3.6, for the same two radial limits listed 
in Table|7] However, these authors have used different distance and 
My for this target. Taking this into account, we have corrected their 
values and obtained Sjv = 3.4 and Sjv = 5.05, somewhat lower 
than those listed in Table|7] 

For NGC 3115, the only p revious estimation of th e cumulative 
specific frequency is that from [Harris & HarrislfcOl ll) . 5'jv ~ 2.0, 
whi ch is very similar to ou r res ults for his galax y. 

iRhode & Zepj ( |2004|) and lHarris & HarrisI l l201lh determine a 
value of Sjv = 1.2 ± 0.3 for NGC 3379, in good agreement with 
our result. 



5 CONCLUSIONS 

Deep GMOS multi-colour photometry has been obtained for the 
GCS of five early type galaxies. The primary results of this study, 
that will be combined with the spectroscopic results for a further 
discussion in a following paper, are: 

• All the studied GCS show bimodal integrated colour- 
distributions. Even the less clear case, NGC 524, appears bimodal 
when analyzed with the RMIX software. The adoption of a mass 
to luminosity ratio {M/L)b= 8 indicates a minimum stellar mass 
close to 1.21O^^M0 for our galaxy sample, i.e., all these systems 
are above the minimu m mass where bim odality becomes detectable 
as a common feature jPeng et al.l2006l) . 

• The mean colour of the peaks of red GCs follow a well de- 
fined luminosity-colour relation in the sense that these peaks be- 
come redder with increasing galaxy luminosity. This is indicative 
of a larger GCs chemical abundance scale for the most luminous 
galaxies which, being more massive, are also able of a more ef- 
ficient enrichment. We note that the most massive system in this 
work (NGC 4649) lies outside the plane defined by GC systems 



in the logarithmic 3-D space determined by galaxy stellar mass, 
projected stellar mass density and GCs formation effic iency, pos- 
sibW indicating a merger past history as suggested by iForte et al] 
ll2009[). The GC blue peaks, in turn, show no significant correla- 
tion or, possibly, only a \ery mild one , with galaxy lum inosity a 
situation already noted in the literature dPeng et al.ll2006l) . 

• An empirical (g'-i') vs [Z/H] relation was obtained 
based solely on GMOS data. A linear fit yields: [Z/H] = 
3.51(±0.21)(5' - i')o - 3.91(±0.20), which provides a good ap- 
proximation, although the eventual presence of a second order term 
at the low metallicity end will require the inclusion of more GCs. 

• A comparison of the colour spreads of both blue and red GCs, 
once transformed to metallicity, clearly indicates that both popula- 
tions are widely different in terms of chemical abundance, a fact 
that when associated with the different spatial distributions of these 
subpopulations, strongly indicates a distinct nature for the two pop- 
ulations. 

• A number of resolved/marginally resolved GC candidates 
were detected in all five galaxies. Some of them are spectro- 
scopically confirmed as GCs in the literature. In principle the 
GCs upper mass sho ws a dependence with galaxy luminosity (e.g. 
IVillegas et all 1201 Ol) and then the detection of resolved clusters 
could be a more frequent situation in these systems. However, fur- 
ther analysis will require the inclusion of distance effects. 

• A population of Ultra Compact Dwarf (UCD) candidates (re- 
solved and unresolved ones with -15 > Mj > -12 ) were also de- 
tected in all the galaxies. NGC 3923 shows the highest number of 
associated objects, followed by NGC4649 and NGC 3115. This re- 
sult indicates that UCDs are indeed a common feature although the 
small size of our sample does not allow a definite conclusion in 
statistical terms. 

• This paper reports the first detection of the so called "blue tilt" 
in NGC 4649 using ground based observations. As a first approach 
we adopted a linear re lation between GC colours and magnitudes 
although, as shown bv lHarrisl ( l2009bh on the basis of a more nu- 
merous GCs sample, the tilt may exhibit some degree of curvature 
that becomes more evident for the brightest clusters. The tilt trans- 
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Table 7. Parameter for fitted Gaussian GCLF. The last four columns give the total number of GC and the Specific Frequency obtained considering limits of 50 
and 100 Kpc for the GCS. 



Galaxy TOM a TOM^^^ ct,,^^ TOMt,„^ a^^e Ntso Sjvbo Ntjoo Sjvioo 



NGC4649 22.95 ± 0.05 1.35 ± 0.05 22.92 ± 0.07 1.28 ± 0.09 22.96 ± 0.07 1.28 ± 0.09 -22.44 3310±I30 3.5±O.I 4690±980 5.0± 1 

NGC 3923 23.84 (fixed) 1.42 ± 0.05 ... ... ... ... -22.03 3019±70 4.1±0.1 4580±820 6.3±l.l 

NGC 3115 21.78 ± 0.08 0.91 ± 0.07 21.59 ± 0.07 0.71 ± 0.05 22.93 ±0.15 1.04 ± 0.12 -21.13 546±80 1.9± 0.2 571±190 2.0±0.7 

NGC 3379 22.46 ± 0.37 0.85 ± 0.23 ... ... ... ... -20.88 216±42 0.9± 0.2 226±100 1.0±0.4 



lates into an approximate mass-metallicity relation, given by Z oc 
jyj-o.28±o.03 xhis result is in excellent a greement withe the mean 
trend obtained bv lCockcroft et al.l ll2009h . 

• All the GCS density profiles are reasonably well fitted by a 
power law, although, in the case of NGC 4649, an R^/* dependence 
produces a better fit in terms of the residuals. As found in other 
galaxies, the red GCs show a more concentrated spatial distribution 
than the blue ones, again lending support to the idea that bimodality 
is a real and distinctive feature and not a mere result of a particular 
colour-metallicity relation. 

• We confirm th e very low spa t ial co ncentration of the NGC 
3923 GCS found bv'Si kkema et al.l l llOOel) . In contrast to these au- 
thors, we find a significant difference of the areal density slopes of 
the red and blue GCs. 

• We have measured radial colour gradients for the blue and 
red GC subpopulations separately, and find evidence for statisti- 
cally significant gradients in several galaxies out to several effective 
radii (tens of kpc). The steeper ones are those present in the two SO 
galaxies included in the sample. These galaxies have smaller stellar 
masses than the brightest ellipticals i n this study showing shallower 
gradients. This is in agreement with lTortora et al. who find 
a non monotonic behaviour of colour gradients with galaxy mass. 
A possible interpretation of this trend assumes that more massive 
galaxies may have experienced mergers that lead to a dilution of 
the chemical gradients. 

• Like in other studied systems, the red GC subpopulation has 
similar colours to the galaxy halo stars in their inner region. This 
is consistent with the idea that most of the galaxy luminosity in 
these regions comes from a diffuse stellar population associated 
with these clusters (e.g., FFG07). 

• The TOM and a were determined for the GCLF of NGC 4649 
and NGC 3115. In the case of NGC 3923, only a was fitted. This 
parameter and the density profiles were used to obtain the total GCs 
populations, and the Specific Frequency Sn of these galaxies. The 
highest ^jv = 4.5 value was obtained for NGC 3923, followed by 
NGC 4649, Sn = 3.6. We have obtained new estimations of the 
NGC 3115 and NGC 3379 specific frequencies; Sn = 2.0 ± 07 
and Sn ~ 1-0 ± 0.4 respectively. 
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Table Al. Photometric catalog of all the GC candidates in the five GCS studied in this work. The magnitudes are from 
psf fitting for the unresolved sources and from aperture photometry for the resolved ones. These magnitudes are extintion 
corrected as indicated in the text. 



Object ID 


RA 


DEC 




Vcad 


a'o 








a f 






notes 


N4649GC2 


12:43:33.74 


11:34:02.5 


1708.009 


53 


.761 


23.8009 


0.023 


23 


.1895 


0.023 


22.8163 


0.024 




N4649GC3 


12:43:33.97 


11:34:01.3 


1685.015 


62 


.249 


22.7809 


0.016 


22 


.0355 


0.008 


21.6613 


0.013 




N4649GC4 


12:43:29.01 


11:34:01.2 


2185.490 


62 


.391 


24.0299 


0.013 


23 


.5035 


0.017 


23.1393 


0.025 




N4649GC5 


12:43:38.04 


11:34:02.4 


1273.952 


55 


.273 


23.7579 


0.028 


22 


.9325 


0.016 


22.5213 


0.019 




N4649GC6 


12:43:41.82 


11:34:02.4 


892.493 


55 


.250 


23.6949 


0.015 


23 


.2035 


0.015 


22.9523 


0.025 




N4649GC7 


12:43:41.24 


11:34:03.2 


950.987 


50 


.004 


23.9169 


0.020 


23 


.3875 


0.019 


23.3253 


0.032 




N4649GC8 


12:43:34.64 


11:33:58.9 


1617.590 


78 


.958 


21.2469 


0.008 


20 


.7345 


0.011 


20.5343 


0.013 




N4649GC9 


12:43:37.67 


11:34:00.1 


1310.816 


70 


.643 


22.5349 


0.010 


22 


.0455 


0.017 


21.7833 


0.022 


a 































The complete photometry for GC candidates is available in the on-line version. 



Table A2. The photometric data for all the UCD candidates in the five GCS studied in this work are given in Table ED 
The magnitudes are from psf fitting for the unresolved sources and from aperture photometry for the resolved ones. These 
magnitudes are extintion coiTected as indicated in the text. 



Object ID 


RA 


DEC 




Vcad 


3(1 


%' 


'■o 


(7 / 


»0 




notes 


N4649UCDI 


12:43:41.23 


11:30:48.3 


952.676 


1388.875 


18.4779 0.003 


17.8995 


0.007 


17.7183 


0.006 







a 



^ The complete photometry for UCD candidate is available in the on-line version. 



